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Over  the  past  two  years  we  have  pursued  actively  new  methodologies  in 
the  mathematical  modeling  and  numerical  simulation  of  the  dynamics  of  flexible 
structures  undergoing  large  overall  motions.  This  two  year  research  effort,  spon¬ 
sored  by  the  AFOSR  under  grant  DJA/AFOSR  86-0292  with  Stanford  Univer¬ 
sity,  has  resulted  in  15  publications,  a  similar  number  of  publications  in  closely 
related  fields  (see  Appendix  II),  and  a  strong  collaboration  (in  a  related 
research  effort  also  sponsored  by  AFOSR)  with  the  groups  of  J.E.  Marsden  at 
Berkeley  (currently  at  Cornell)  and  P.S.  Krishnaprasad  at  Maryland^. 


First,  we  outline  in  Section 
2  the  scope  of  the  proposed  research,  and  summarize  the  research  goals  as 
stated  in  the  original  proposal.  In  Section  3  we  provide  a  brief  summary  of  the 
objectives  accomplished  during  the  past  two  year  effort.  Explicit  reference  is 
made,  were  appropriate,  to  the  work  published  under  the  present  research  con¬ 
tract.  We  remark  that  all  of  the  proposed  objectives  (summarized  in  Section  2) 
have  been  accomplished.  Section  4  is  concerned  with  additional  topics  pursued 
in  these  past  two  years,  which  which  fall  within  the  subject  area  of  this  research 
and  are  not  explicitly  identified  in  the  original  proposal.  The  contents  of  these 
three  sections  is  essentially  a  summary  of  the  results  reported  in  the  final  report 
submitted  to  AFOSR  in  spring  of  1987.  We  refer  to  this  latter  document  for 
further  details  and  elaboration  on  the  nature  of  these  results. 


2.  Review  of  scope  and  objectives  of  the  original  AFOSR  proposal 

As  stated  in  the  abstract  of  our  original  AFOSR  proposal,  the  ongoing 
research  under  the  present  AFOSR  contract  constitutes  an  effort  towards  the 
development  of  modeling  and  computational  tools  for  the  integrated  design  of 
flexible  structures  undergoing  large  overall  motions  Approaches  in  wide  use  at 
present  rely  crucially  on  the  use  of  a  floating  reference  frame  that  follows  the 
overall  rigid  body  motion  of  the  system  thus  allowing  the  use  of  simplified 
structural  models  based  on  the  assumption  of  "small  strains."  From  a  physical 
standpoint  this  approach  precludes  the  modeling  of  physically  important 


t  An  example  of  this  collaboration  is  the  forthcoming  AMS  conference  in  Maine,  joint- 

ly  organized  by  P.S.  Krishnaprasad,  J.E.  Marsden  &  J.C.  Simo.  J 
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nonlinear  effects,  such  as  large  deformations,  nonlinear  vibrations  and  struc¬ 
tural  damage  that  may  arise,  for  instance,  as  a  result  of  control  system  break¬ 
down.  Thus,  worse  case  situations  cannot  be  accounted  for  in  the  design  pro¬ 
cess.  From  a  computational  stand  point,  approaches  currently  in  use  lead  to 
discrete  equations  of  motion  that  are  highly  coupled  in  the  inertia  terms  as  a  result 
of  inertia  effects  due  to  rotation  of  the  floating  frame. 

2.1.  Overview  and  general  remarks 

The  central  theme  of  the  present  work  is  the  development  of  an  a  new 
alternative  methodology  that  hinges  on  the  use  of  geometrically  exact  structural 
models  to  mathematically  characterize  the  flexibility  properties  of  the  structure. 
This  approach  enjoys  the  following  advantages  over  traditional  formulations 

i.  Geometrically  exact  models  possess  full  50(3)  invariance;  hence,  they  can 
be  used  directly  in  the  modeling  of  coupled  rigid  body-flexible  appendage 
systems  without  resorting  to  the  introduction  of  the  so-called  floating 
frames.  This  is  of  particular  importance  in  applications  concerned  with  the 
dynamics  of  orbiting  satellites  with  flexible  appendages.  In  fact,  the  entire 
dynamics  of  the  coupled  system  can  be  referred  directly  to  the  inertial 
frame. 

iL  One  has  the  added  capability  of  modeling  important  physical  phenomena 
precluded  by  linearized  theories;  such  as  large  deformation,  bifurcation  and 
loss  of  stability.  In  addition,  complete  inertia  effects  can  be  exactly 
accounted  for  without  simplification.  Such  a  generality  is  warranted  in 
many  engineering  applications;  for  instance,  helicopter  and  rotor  blades, 
modern  highly  flexible  robot  arm  manipulators,  and  space  structures. 
Finally,  nonlinear  geometric  effects,  which  play  an  important  role  in  the 
analysis  of  rotating  structures,  are  automatically  accounted  for.  In  this 
regard,  it  should  be  noted  that  use  of  linearized  theories  in  the  dynamic 
analysis  of  rotating  structures  may  lead  to  completely  erroneous,  non¬ 
physical,  results. 

iu.  Computationally,  by  referring  the  motion  of  the  system  to  the  inertial 
frame,  extensive  coupling  in  the  inertia  operator  due  to  Coriolis,  centrifu¬ 
gal,  and  inertia  effects  due  to  rotation  of  the  floating  frame  is  avoided.  A 
case  in  point  concerns  the  dynamics  of  planar  systems  for  which  the  inertia 
operator  becomes  linear.  Furthermore,  because  of  the  full  S0( 3)- 
invariance  property,  the  proposed  formulation  is  specially  well-suited  for 
the  numerical  simulation  of  multi-body  interconnected  systems  in  complex 
configurations.  This  is  the  case,  for  instance,  of  closed-loop  chains  under¬ 
going  large  overall  motions. 
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2.2.  Objectives  of  the  investigation 

The  computational  and  physical  considerations  briefly  summarized  above 
led  to  the  following  research  goals  to  be  accomplished  in  the  course  of  a  two- 
year  investigation  i 

1.  Development  of  a  new  class  of  covariant  time  stepping  algorithms  for  the 
numerical  simulation  of  the  transient  response  of  fully  nonlinear  flexible 
structures.  For  this  purpose,  concepts  of  modern  differential  geometry 
play  a  central  role. 

2.  Application  of  the  methodology  proposed  above  to  the  dynamics  of  multi¬ 
body  flexible  systems. 

3.  Development  of  a  new  class  of  fully  nonlinear  formulations  for  the  dynam¬ 
ics  of  structural  elements  such  as  plates  and  shells. 

4.  Extension  to  more  general  constitutive  behavior  that  includes  possibly 

nonlinear  dissipative  mechanisms,  and  structural  damage.  Formulation  of 
general  integration  algorithms  for  numerical  simulation  of  this  constitutive 
behavior.  I 

5.  Development  of  effective  numerical  algorithms  for  the  simulation  of  orbit¬ 
ing  LSS  composed  of  beam-like  elements.  An  essential  ingredient  is  the 
separation  and  effective  treatment  of  the  far  field  and  the  near  field. 

6.  Simulation  and  design  of  flexible  structures  undergoing  large  overall 
motions.  This  includes  fast  slewing  and  precision  pointing  in  a  wide  range 
of  applications,  and  appropriate  treatment  of  follower  and  non-conservative 
loading  (e.g.  gravity,  actuator  control  forces  on  large  space  structures). 

A  summary  of  the  work  completed  during  the  past  two  year  effort,  and  its  rela¬ 
tion  to  the  aforementioned  research  goals  is  given  in  the  following  section. 

3.  Summary  of  objectives  proposed  and  accomplished  in  the  two  years 

All  of  the  tasks  proposed  in  the  original  proposal  to  the  AFOSR  outlined 
above  have  been  completed  in  the  two-year  period  spanned  by  our  investiga¬ 
tion.  An  outline  of  our  results  follows. 

3.1.  Time-stepping  covariant  algorithms  for  nonlinear  structural  dynamics. 

The  methodology  for  planar  systems  developed  in  Simo  &  Vu-Quoc 
[1986a,b],  is  first  extended,  in  Simo  &  Vu-Quoc  [1986c],  within  the  context  of 
the  static  problem,  to  the  full  three-dimensional  case.  The  general,  fully  cou¬ 
pled  dynamic  case  is  considered  in  Simo  &  Vu-Quoc  [1988a].  A  summary  of 
the  methodology  is  given  in  the  invited  contribution  in  Simo  &  Vu-Quoc 
[1987].  Some  salient  features  of  this  work  are  the  following: 
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i.  A  complete  convergence,  and  accuracy  analysis  is  given. 

iL  Methods  of  differential  geometry  are  rigorously  used  to  exploit  the  under- 
lying  geometric  structure  of  the  model.  In  particular,  configurations  are 
updated  by  means  of  the  discrete  version  of  the  exponential  map  in 
R3xS0(3),  whereas  velocity  and  accelerations  are  updated  by  making  use 
of  the  discrete  parallel  transport  in  R3xS0(3). 

iti.  The  algorithmic  procedure  can  accommodate  any  degree  of  geometric  non* 
linearity.  In  particular,  limit  and  bifurcation  points  are  traverse  by  con¬ 
tinuation  procedures,  as  discussed  in  Simo,  Wriggers,  Schweizerhof  &  Tay¬ 
lor  [1986]. 

iv.  The  formulation  and  numerical  implementation  includes  actuator  and  fol¬ 
lower  forces  of  importance,  for  instance,  in  attitude  control  by  means  of 
thrusters. 

v.  A  singularity  free  parametrizations  of  the  rotation  field  based  on  the  use  of 
quaternion  parameters  is  employed,  which  appears  to  be  optimal  for  large- 
scale  computations. 

A  numerical  simulation  that  illustrates  the  performance  of  the  methodol¬ 
ogy  described  above  is  contained  in  Figure  1. 


Figure  1.  Rod  undergoing  large  amplitude  vibration  superposed  onto  three  di¬ 
mensional  overall  motions. 


3.2.  Application  to  multibody-dynamics. 

The  numerical  simulations  presented  in  Simo  &  Vuoc  [1987,1988a],  and  in 
Vu-Quoc  and  Simo  [1987]  demonstrate  the  applicability  of  the  proposed  metho¬ 
dology  to  multibody  dynamics.  As  an  illustration.  Figure  2  shows  the  simula¬ 
tion  of  the  dynamics  of  an  all-flexible  closed  loop  chain  system  undergoing  large 
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overall  motions;  an  example  of  a  class  of  problems  to  which  the  highest  degree 
of  difficulty  is  assigned  in  the  well-known  classification  of  Ho  &  Herber  [1985]. 
Furthermore,  the  example  exhibits  an  added  degree  of  difficulty  which  appears 
to  have  no  counterpart  in  the  current  literature  on  multibody  dynamics:  no  res - 
triction  is  placed  on  the  degree  of  flexibility  of  the  chain  members. 


Figure  2.  AU-flexiblc,  dosed  loop  chain  undergoing  finite  vibration  and  large 
overall  motions. 


The  fact  that  classical  methodologies  relying  on  ad-hoc  linearized  models 
can  result  in  completely  erroneous  non-physical  results  is  demonstrated,  both 
analytically  and  numerically,  in  Simo  &  Vu-Quoc  [1988b].  This  analysis  also 
provides  a  systematic  means  of  developing  higher  order  theories  from  the 
geometrically  exact  theory  by  an  asymptotic  expansion  procedure. 

3.3.  Well-conditioned  methodologies  for  satellite  dynamics 

That  flexibility  plays  a  significant  role  in  the  stability  of  satellites  is  a  well- 
established  fact  that  was  recognized  as  early  as  the  start  of  the  space  program,  f 
In  Vu-Quoc  &  Simo  [1987],  by  exploiting  the  full  SO  (3) -invariance  of  the  pro¬ 
posed  geometrically  exact  models,  the  dynamics  of  the  orbiting  satellite  is 
referred  directly  to  the  inertial  frame.  This  point  of  view  represents  a  radical 
departure  from  traditional  approaches,  and  offers  the  following  advantages  with 
respect  to  traditional  formulations  which  rely,  typically,  on  floating  frame  sys¬ 
tems  and  make  crucial  use  of  the  "small  deformations”  assumption  (e.g.,  Cana- 
vin  &  Likins  [1977]), 


t  As  noted  in  the  original  proposal  submitted  to  AFOSR,  the  instability  of  Explorer  I, 
the  first  U.S.  satellite,  was  caused  by  the  flexibility  of  its  small  wire  turnstile  antennae. 
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i.  From  a  computational  point  of  view,  the  task  of  coupling  rigid  body  and 
flexible  appendage  becomes  remarkably  simple. 

iL  The  structure  of  the  inertia  terms  simplifies  considerably. 

iu.  Large  deformations,  which  are  expected  to  take  place  in  very  flexible 
appendages  and  are  associated  with  important  physical  phenomena  such  as 
loss  of  stability  and  buckling,  can  be  predicted.  , 

*  '  i 

UL  Numerical  ill-conditioning  associated  with  the  enormous  differences  in  ord¬ 
ers  of  magnitude  between  orbital  radius  and  amplitude  of  deformation  is 
avoided  by  referring  the  dynamics  of  the  orbiting  satellite  to  a  parallel 
translate  of  the  inertial  frame  located  in  the  vicinity  of  the  coupled  system 
(not  necessarily  at  the  center  of  mass  of  the  system). 

iv.  In  addition  to  avoiding  ill-conditioning,  the  aforementioned  technique  has 
the  advantage  of  effectively  decoupling  the  dynamics  of  the  system  into  a  far 
field  and  a  near  field  problems. 

v.  We  remark  that  the  formulation  can  accommodate  configuration  dependent 
follower  loads.  This  feature  is  of  crucial  practical  importance  in  attitude 
control  of  orbiting  satellites  since  "gravity  gradient"  and  actuator  or  control 
forces;  such  as  gas  jets  or  ion  thrusters,  follow  the  deformation  of  the 
space  structure 

3.4.  General  nonlinear  dynamic  rod  models 

A  main  limitation  of  existing  geometrically  exact  mathematical  models  for 
beam-like  structures,  sec  Le.,  Amman  [1974],  Antman  &  Jordan  [1975],  Reiss- 
ner  [1973,1981],  is  their  inability  to  accommodate  the  effect  of  torsional  warping 
distortion  of  the  cross  section  of  the  rod.  This  effect  plays  a  dominant  role  in 
many  situations  of  engineering  interest;  particularly  those  concerned  with  light¬ 
weight  structures  that  employ  thin-walled  elements  with  open  or  closed  channel 
sections.  A  case  in  point  is  the  design  of  robot-arm  manipulators  where  light¬ 
weight  consideration  preclude,  almost  always,  the  use  of  a  solid  cross-section. 
Partly  because  of  these  considerations,  the  problem  of  coupling  bending-torsion 
has  received  considerable  attention  in  recent  years;  Le.,  Reissner 
[1979, 1983a,  1983b,  1984],  Hodges  [1980],  Krenk  [1983a, b],  Krenk,  &  Gun- 
neskov,  [1985].  No  results,  however,  appear  to  exist  in  a  geometrically  exact 
context. 

In  Simo  &  Vu-Quoc  [1988c],  we  have  extended  the  model  proposed  in 
Simo  [1985]  to  accommodate  the  effect  torsional  warping  distortions  of  the  cross 
section.  In  addition,  our  earlier  computational  framework  has  been  extended  to 
accommodate  the  effect  of  warping.  A  numerical  simulation  that  illustrates  the 
important  effect  of  warping,  and  exhibits  the  full-nonlinear  capabilities  of  the 
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Figure  3.  Lateral  loss  of  stability  and  bifurcation  of  a  right-angle  frame  is 
presented  illustrating  the  effect  of  coupled  torsion-bending-warping  in  the  bifurca¬ 
tion  diagram. 


3.5.  Formulation  of  geometrically  exact  plate  and  shell  models 

Over  the  past  two  decades  engineering  shell  analysis  has  been  dominated, 
to  a  large  extent,  by  the  so-called  degenerated  solid  approach,  which  finds  its 
point  of  departure  in  the  paper  of  Ahmad,  Irons  &  Zienkiewicz  [1970].  The 
works  of  Ramm  [1977],  Parish  [1981],  Hughes  St  Liu  [1981a,b],  Hughes  St 
Camoy  [1983],  Bathe  &  Dvorkin  [1984],  Hallquist,  Benson  St  Goudreau 
[1986],  Parks  St  Stanley  [1986],  and  Liu,  Belytschko,  Law  St  Lam  [1987], 
among  many' others,  constitute  representative  examples  of  this  methodology. 
Our  work  in  Simo  St  Fox  [1988],  on  the  other  hand,  constitutes  a  radical  depar¬ 
ture  from  the  aforementioned  methodology.  In  a  sense,  the  proposed  approach 
represents  a  return  to  the  origins  of  classical  nonlinear  shell  theory,  which  has 
its  modern  point  of  departure  in  the  fundamental  work  of  the  Cosserats  [1909], 
subsequently  rediscovered  by  Ericksen  St  Truesdell  [1958],  and  further  ela¬ 
borated  upon  by  a  number  of  authors;  notably  Green  St  Laws  [1966],  Green  & 
Zerna  [I960],  Cohen  &  DaSilva  [1966],  Naghdi  [1972,1980],  Green  &  Naghdi 
[1974],  and  Antman  [1976a,b]. 

In  Simo  &Fox  [1987],  we  have  initiated  the  mathematical  formulation  and 
optimal  parametrization  of  a  general  class  of  nonlinear,  geometrically  exact 
plate  and  shell  models,  which  is  amenable  to  numerical  analysis  and  suitable  for 
large-scale  computation.  In  Simo,  Marsden  St  Krishnaprasad  [1987]  we  have 
undertaken  what  appeals  to  constitute  the  first  mathematical  analysis  of  the 
Hamiltonian  structure  underlying  this  important  class  of  models.  In  contrast 


J.C.  Sitno 


10 


with  classical  treatments,  the  emphasis  of  this  work  is  on  the  underlying 
geometric  structure  of  the  nonlinear  shell  theory.  The  understanding  and 
exploitation  of  this  structure  leads  to  the  following  noteworthy  results 


Figure  4.  Large  deformation  and  torsional  collapse  of  a  cylindrical  shell.  Actu¬ 
al  deformations,  shown  without  magnification,  are  attained  in  Ihrce  load  steps. 

i.  Optimal  singularity-free  parametrizations  of  the  director  (rotation)  field  arc 
developed  by  identifying  the  geometric  structure  of  the  configuration  space 
C;  essentially  a  differentiable  manifold  of  mappings  modeled  on  the  pro¬ 
duct  manifold  S 2  x  1R3,  where  S 2  denotes  the  unit  sphere. 

u.  Update  procedures  for  the  director  and  associated  rotation  field,  which  are 
exact  and  preserve  objectivity  for  any  magnitude  of  the  director  and  rota¬ 
tion  increment,  are  constructed  by  exploiting  the  following  two  geometric 
objects:  (a)  A  discrete  version  of  the  exponential  map  in  the  unit  sphere 
S2,  and  (b)  The  realization  of  the  unit  sphere  as  the  non-trivial  sundle 
SO(3)/SK  This  realization  corresponds  to  a  reduction  of  SO  (2)  by  the 
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group  of  rotations  in  the  plane  S*. 

ill.  Exact  enforcement  of  constraints  such  as  in  extensibility  of  the  director  field  is 
achieved  exactly  without  resorting  to  ad-hoc  procedures  (or  penalty 
methods)  by  exploiting  the  aforementioned  realization  of  the  unit  sphere 
as  the  bundle  SO(3)/Sl. 

iv.  A  parameterization  of  the  basic  weak  form  of  the  momentum  equations  is 
given  which  avoids  explicit  appearance  of  the  Riemannian  connection  on 
the  mid-surface.  Objects  such  as  Christoffel  symbols  and  covariant  deriva¬ 
tives,  which  are  not  readily  accessible  in  a  computational  context,  never 
appear  explicitly  in  our  formulation.  Consequently,  a  major  drawback  in 
the  numerical  treatment  of  of  classical  shell  theory  is  entirely  bypassed. 

Some  preliminary  numerical  simulations  which  illustrate  the  effectiveness 
of  our  methodology  are  given  in  Figure  4  and  Figure  5. 


Figure  5.  Large  deflection  of  a  quarter  of  a  cylindrical  shell.  The  deformation 
shown  (with  no  magnification)  is  obtained  in  one  single-step  bading. 


3.6.  Numerical  analysis  of  geometrically  exact  shell  models 

A  point  frequently  made  concerning  the  degenerated  solid  approach  is  that 
it  avoids  the  mathematical  complexities  associated  with  classical  shell  theory, 
and  hence  is  better  conditioned  for  numerical  implementation.  In  Sinio,  Fox  & 
Rifai  1 1988a, bj,  on  the  other  hand,  we  have  demonstrated  that  stress-resultant 
geometrically  exact  shell  models,  phrased  in  the  geometric  setting  discussed  in 
Section  3.5  are  indeed  ideally  suited  for  large  scale  numerical  simulations,  and 
amenable  to  numerical  analysis.  In  order  to  evaluate  the  proposed  finite 
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element  formulation  and  assess  its  performance  in  standard  benchmark  prob¬ 
lems,  we  have  initially  restricted  our  attention  to  the  linear  theory  which  is 
obtained  by  consistent  linearization  of  the  full  nonlinear  formulation  in  Simo  Sc 
Fox  [1988].  This  is  a  preliminary  step  which  must  be  taken  before  addressing 
the  numerical  analysis  aspects  involved  in  the  implementation  of  the  full  non¬ 
linear  theory.  Some  of  the  results  emanating  from  our  recent  work  are  as  fol¬ 
lows  • 


Figure  6.  The  pinched  hemisphere:  A  standard  benchmark  problem  for  inexten - 
sional  response  and  large  rigid  body  morions. 


L  The  underlying  geometric  structure  of  the  continuum  theory  is  preserved 
exactly  by  our  discrete  finite  element  interpolation  of  the  director  and  the 
linearized  director  fields.  A  discrete  singularity-free  mapping  between  the 
five  and  the  six  degree  of  freedom  formulation,  which  appears  to  be 
canonical,  is  constructed  again  by  exploiting  the  realization  of  S2  as  the 
bundle  50(3)/S1.  As  a  result,  :he  proposed  numerical  formulation 
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exhibits  excellent  performance  under  large  rigid  body  motions.  This  is 
illustrated  in  Figure  6  in  a  classical  benchmark  problem.  In  fact,  our  bil¬ 
inear  finite  element  formulation  exhibits  performance  superior  to  current 
higher  order  state-of-the-art  finite  element  formulations  based  on  the 
degenerated  solid  concept. 


Figure  7.  Convergence  results  for  the  rhombic  plate:  A  severe  test  of  the  perfor¬ 
mance  of  plate  bending  elements. 


The  proposed  mixed  interpolation  for  the  bending  field,  which  emanates 
from  a  Hellinger-Reissner  variational  formulation,  leads  to  a  finite  element 
formulation  which  is  free  from  spurious  zero  energy  modes,  and  appears  to 
yield  optimal  results  in  standard  benchmark  problems,  (Le.,  MacNeal  & 
Harder  [1984]).  For  the  classical  rhombic  plate  test  (a  singularity  dom¬ 
inated  problem)  for  instance  (see  Figure  7)  the  performance  of  the  ele¬ 
ment  is  clearly  superior  to  existing  four-node  and  even  nine-node  ele¬ 
ments;  including  T1  (Hughes  &  Tezduyar  [1981]),  the  equivalent  element 
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of  Bathe  St  Dvorkin  [1984],  and  the  four  and  nine-node  ANS  formulations 
of  Parks  St  Stanley  [1986]. 


Figure  8.  Convergence  results  for  Cook’s  membrane  problem:  A  classical 
benchmark  for  membrane  dominated  problems . 


id.  The  proposed  interpolation  for  the  membrane  field,  also  based  on  a  mixed 
Hellinger-Reissner  formulation,  appears  to  yield  optimal  results  in  standard 
benchmark  tests  for  membrane  dominated  problems;  such  as  Cook’s  mem¬ 
brane  test,  see  Figure  8. 

rv.  Although  our  finite  element  interpolations  satisfy  exactly  all  the  standard 
engineering  test  for  convergence;  in  particular  the  modified  version  of  the 
patch  test  (Taylor,  Simo,  Zienkiewicz  St  Chan  [1986])  no  rigorous 
mathematical  proof  of  convergence  is  yet  available.  Nevertheless,  we  note 
th&i  with  the  exception  of  linear  analysis  of  Koiter’s  model  by  Bernadou  St 
Boisserie  [1982],  and  the  recent  work  of  Brezzi  St  Fortin  [1985]  and 
Arnold  St  Falk  [1987]  for  the  plate  problem ,  no  rigorous  convergence 
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results  for  the  discrete  shell  problem  appear  to  be  available  in  the  numeri¬ 
cal  analysis  literature. 

4.  Additional  research  accomplished  in  the  two  year  effort 

In  addition  to  the  research  objectives  outlined  in  the  preceding  section,  two 
complementary  lines  of  investigation  have  been  explored  the  results  of  which 
are  summarized  below.  The  theoretical  results  obtained  in  the  first  area  have 
proven  essential  in  the  formulation  and  numerical  analysis  aspects  of  the 
geometrically  exact  shell  models  described  above. 

4.1.  Hamiltonian  structure  of  elasticity,  rods  and  plates 

It  is  now  widely  accepted  that  a  thorough  understanding  of  the  mathemati¬ 
cal  underpinnings  of  elasticity  is  crucial  to  its  analytical  and  computational 
implementation.  In  particular,  in  the  work  reported  above,  a  geometric  under¬ 
standing  of  elasticity  has  proven  essential  in  the  formulation  of  efficient  algo¬ 
rithms  for  geometrically  exact  rod  models.  Underlying  this  geometric  view  of 
elasticity  there  exists  a  remarkably  simple  non-canonical  Hamiltonian  structure 
which  is  explored  and  developed  in  detail  in  Simo,  Marsden  &  Krishnaprasad 
[1987].  Basic  motivations  for  pursuing  this  line  of  investigation  are  the  follow¬ 
ing 

1.  From  a  fundamental  standpoint,  a  geometric  formulation  of  the  equations 
of  elasticity  in  Hamiltonian  form  is  extremely  useful  in  the  analysis  of  a 
wide  range  problems  of  physical  interest.  Typical  examples  include  the 
study  of  nonlinear  stability  of  solids,  fluids,  and  plasmas;  see  Le.,  the 
comprehensive  survey  of  Holm,  Marsden,  Ratiu  .&  Weinstein  [1986];  and 
the  analysis  of  bifurcation  phenomena;  i.e.,  Golubitsky  &  Stewart  [1986], 
and  Marsden  &  Ratiu  [1986]. 

2.  From  a  computational  point  of  view,  geometric  methods  play  an  important 
role  in  the  the  development  of  numerical  schemes  in  nonlinear  elasticity. 
For  instance,  as  already  noted,  the  algorithmic  counterparts  of  the 
exponential  map  and  the  parallel  transport  are  the  essential  tools  in  the  for¬ 
mulation  of  algorithms  discussed  in  Section  3.1.  Similar  techniques  are 
useful  in  nonlinear  inelasticity;  Le.,  Simo  [1987a,b,c]. 

3.  A  geometric  approach  along  with  the  use  of  Hamiltonian  formalism  has 
proven  useful  in  successfully  tackling  problems  involving  coupled  rigid 
bodies/ flexible  attachments.  An  example  of  this  is  the  rigorous  stability 
analysis  of  Krishnaprasad  &  Marsden  [1987]  for  i  restricted  class  of  flexi¬ 
ble  appendages  coupled  with  a  rigid  body. 

With  this  motivation  at  hand,  the  results  obtained  in  Simo,  Marsden  & 
Krishnaprasad  [1988]  may  be  briefly  summarized  as  follows 
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L  A  complete  and  systematic  development  of  the  Hamiltonian  structure  of 
nonlinear  elasticity,  and  geometrically  exact  rods  and  shells  models  is 
given.  An  interesting  point  of  this  development  is  the  central  role  played 
by  the  convective  representation  which,  for  rigid  bodies,  collapses  to  the 
familiar  description  of  the  motion  in  body  coordinates. 

iL  Identification  of  an  underlying  non-canonical  Lie-Poisson  structure  com¬ 
mon  to  all  geometrically  exact  models,  and  emanating  from  the  convective 
representation  is  given.  A  complete  derivation  of  the  corresponding  Lie- 
Poisson  brackets  is  also  given.  These  Lie-Poisson  structures  will  become 
essential  in  future  stability  calculations  of  coupled  rigid  bodied  with  flexible 
appendages.  A  systematic  analysis  of  the  dynamics  of  geometrically  exact 
models  is  already  in  progress  and  the  results  are  reported  in  Simo, 
Marsden  &  Posbergh  [1988]. 

nr.  For  the  case  of  geometrically  exact  shell  models,  a  novel  geometric  picture 
emanates  from  this  analysis  which  is  currently  perceived  as  playing  a  cen¬ 
tral  role  in  future  algorithmic  treatments.  In  particular,  as  the  discussion 
of  Section  3.5  and  Section  5.6  suggests,  the  evolution  of  the  director  field 
in  a  one-director  Cosserat  surface  is  intimately  connected  to  the  geometric 
structure  associated  with  the  group  of  mappings  from  an  open  set  QcR2 
onto  the  unit  S2-sphere.  ! 

I 

4.2.  Improved  material  modeling.  Material  damage  and  inelasticity 

As  pointed  out  in  our  AFOSR  proposal,  future  generations  of  space  struc¬ 
tures  are  expected  to  exhibit  highly  nonlinear  constitutive  response.  Weight 
and  strength  considerations  motivate  a  growing  use  of  composite  materials, 
such  as. the  traditional  graphite/ epoxi  composites,  or  the  relatively  recent  tech¬ 
nology  of  metal-matrix  composites.  Effective  use  of  this  latter  class  of  materi¬ 
als  often  requires  the  matrix  to  operate  in  a  semi-plastic  regime,  leading  to  a 
highly  nonlinear  response.  By  contrast,  the  nonlinearity  in  the  response  of 
graph ite-epoxi  composites  emanates  from  a  physically  different  mechanism 
related  to  damage  degradation  of  the  elastic  properties  of  the  materiai 

Because  of  the  complexity  of  the  physical  processes  involved,  which  take 
place  at  the  micro-structural  level,  characterization  of  damage  degradation  and 
plastic  flow  is  often  made  by  means  of  phenomenolo|gical  constitutive  models. 
This  has  lead  to  the  expanding  field  of  damage  mechanics.  In  its  present  stage, 
however,  a  solid  mathematical  foundation  appears  appears  to  be  lacking. 

Our  research  efforts  in  this  direction  aim  at  the  development  of  mathemat¬ 
ically  tractable  constitutive  models  for  inelasticity  and  damage  response  which 
are  suited  for  large  scale  computation.  The  work  in  Simo  [1987a,b,c],  and 
Simo  &  Ju  [1987a,b,c]  constitutes  an  effort  in  this  direction  for  plastic. 
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viscoelastic-damage,  and  elastic-damage  response.  The  work  in  Simo  &  Honein 
[1988],  constitutes  a  systematic  development  of  a  yariational  formulation  of 
elastoplasticity  along  with  conservation  laws  within  the  context  of  Noether’s 
theorems.  Relevant  numerical  analysis  and  computational  aspects  are  addressed 
in  Simo  &  Taylor  [1986],  Simo  &  Hughes  [1986,1987],  Simo  Kennedy  & 
Govindjee  [1988],  Simo,  Kennedy  &  Taylor  [1988],  and  Simo  [1987d]. 


Appendix  I 

Presentations  at  Meetings  and  Conferences  in  1986-1988. 

1.  IFAC  conference,  July  1986. 

2.  World  congress  on  Computational  Mechanics,  October  1986. 

3.  ASME  Winter  Anual  Meeting,  December,  1986. ; 

4.  International  Conference  in  Numerical  Methods  Theory  and  Applications 
(NUMETA),  July,  1987 

5.  Workshop  on  Dynamics  and  Controls,  (Sponsored  by  NASA-SDIO)  Jet 
Propulsion  Laboratory ,  September  1987. 

6.  IFAC  conference,  December  1987. 

7.  International  Conference  on  Computational  Engineering  Science,  April, 
1988.  Key-Note  address. 

8.  Instituto  per  le  Applicazioni  del  Calcolo,  Rome,  Italy,  May  1988.  (2) 
Invited  Lectures. 

9.  Work-Shop  on  Symmetry  Groups  and  Invariance  in  Continuum  Mechan¬ 
ics.  Cornell  University,  June  1988.  Key-Note  address. 
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10.  International  Meeting  on  Engineering  Science,  Berkeley,  June  1988. 

11.  Summer  Conference  on  Dynamics  and  Multibody  Systems,  August,  1988, 
Organized  by  P.S.  Krishnaprasad,  J.  Marsden  and  J.C.  Simo.  Sponsored  by 
AMS  and  SIAM.  Key-Note  address. 

12.  IFAC  Conference,  December,  1988. 
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